Draft version February 1, 2008 

Preprint typeset using I£T^X style cmulateapj v. 10/09/06 



o 
o 

bJQ' 

< 

00 



SPECTRAL MAPPING RECONSTRUCTION OF EXTENDED SOURCES 

J.D.T. Smith 1 , L. Armus 2 , D.A. Dale 3 , H. Roussel 4 , K. Sheth 2 , B.A. Buckalew 5 , T.H. Jarrett 5 , G. Helou 5 , and 

R.C. Kennicutt, Jr. 6,1 

Draft version February 1, 2008 

ABSTRACT 

Three dimensional spectroscopy of extended sources is typically performed with dedicated integral 
field spectrographs. We describe a method of reconstructing full spectral cubes, with two spatial 
and one spectral dimension, from rastered spectral mapping observations employing a single slit in a 
traditional slit spectrograph. When the background and image characteristics are stable, as is often 
achieved in space, the use of traditional long slits for integral field spectroscopy can substantially 
reduce instrument complexity over dedicated integral field designs, without loss of mapping efficiency 
— particularly compelling when a long slit mode for single unresolved source followup is separately 
required. We detail a custom flux-conserving cube reconstruction algorithm, discuss issues of extended 
source flux calibration, and describe Cubism, a tool which implements these methods for spectral maps 
obtained with ther Spitzer Space Telescope's Infrared Spectrograph. 

Subject headings: methods: data analysis — techniques: spectroscopic — infrared: general 
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1. INTRODUCTION 

Spectroscopy has long been at the forefront of as- 
tronomical research. From uncovering the elemental 
abundances underlying the solar spectrum, to modern 
massively parallel surveys probing the scale and struc- 
ture of the universe, the dispersion of light affords 
an almost preternatural power of discovery at a dis- 
tance. The recent advent of large format optical and 
infrared detector arrays has revolutionized spectroscopy, 
enabling many novel methods which capitalize on the 
growing pixel count. Cross-dispersion of the dispersed 
beam enables higher spectral resolution and longer 



slits i n compact instruments (e.g . Allington- Smit h et al 
TOSI iTull et al.llT9 95: Smi th et alJ Il998t iMcLean et al 



19981 : IWilson et all 12001( 1 The use of multiple slits 
or fibers makes it possible to obtain spectra for hun- 
dreds of objects simult aneously oyer wide fields (e.g. 
Heacox fc Colmesl Il99l lYork et all I2000T: iColless etatl 
20011 iLe Fevre et all 12001b iHook et al.l l2003f ) . Integral 
field spectroscopy of moderately wide areas using image 
slicing, densely packed fiber bundles, o r micro- lens arrays 
efficiently map extended objects (e.g . 
Contend H991 iMandel etaD llOOOt 



Tecza et alJfl 998: 
Dubbcldam et all 



2000t iBacon et al.l 120011 : iPoglitsch et all l2006fl . These 
methods, along with Fabry-Perot scanning and other 
imaging interferometers, fall under the broad category 
of three dimensional spectroscopy, which encompasses a 
set of techniques and instrumental designs which enable 
the observer to build up spectral cubes with two spatial 
and one spectral dimension. 

Given finite detector area, a tradeoff among field of 
view, spectral resolution, and instantaneous wavelength 
coverage must be made for a given three dimensional 
spectrograph design. In the limit of zero background, 
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identical throughput and detector performance, and an 
unchanging point spread function (PSF), any combina- 
tion of these three which utilizes the same fractional 
detector area will offer identical spectral mapping effi- 
ciency. Modern integral field spectrographs typically fa- 
vor larger fields at the cost of resolution or wavelength 
coverage, but have the notable advantage of capturing 
the entire field at once, which mitigates systematics due 
to variations in sky background and transmission, and 
time var iations in the PSF (e.g. in adaptive optics appli- 
cations, iThatte et ai1l2000D . For observing single unre- 
solved sources, however, integral field spectrographs are 
relatively inefficient, since they dedicate a large fraction 
of their detector budget to unused background area. For 
isolated point sources, single slit spectroscopy still main- 
tains an advantage. This advantage is even more acute 
in the ground-based infrared, where, due to the large and 
variable sky brightness, source and sky must be switched 
rapidly throughout the observation — on timescales of 
minutes at near infrared (NIR) wavelengths (1-2.5 (im) 
and up to 10Hz at thermal mid-infrared (MIR) wave- 
lengths (4-20 fim) . In this regime, single long slits of- 
fer another compelling advantage: the ability to chop 
the source between two positions along the slit, perform- 
ing in situ sky subtraction with well-characterized slit 
throughput, and minimal loss of observing efficiency. 

Unfortunately, single long slits are of more limited use 
for three dimensional spectroscopy of extended sources 
from the ground. In particular, variations in the sky 
brightness and transmission, as well as the PSF dur- 
ing the (potentially long) series of integrations required 
to map out large extended regions pose significant chal- 
lenges, especially for sources with surface brightness com- 
parable to or below that of the sky. Despite these draw- 
backs, the scanning method has been used at optical 
and near-infrared wavelengths t o construct moderate- 
sized maps of bright sources (e.g. iBurton fc Afie^ll99l : 
Monteiro et al.l I2005D. and in solar astronomy (e.g. 



Johanneson et afl lT992). but has been eclipsed in recent 



years by dedicated integral field techniques. Though they 
offer high mapping efficiency for extended sources, the 
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reduced efficiency of integral field spectrographs for un- 
resolved sources often dictates that additional long-slit 
modes are included in instrument which offer these ca- 
pabilities, increasing complexity. 

In space, the various limitations of long slits for con- 
ducting three dimensional spectroscopy over wide fields 
are minimized, since the astrophysical background is typ- 
ically smooth and unvarying, and the imag e characteris- 
tics fi xed. The Infrared Spectrograph (IRS, Houck ct al. 
2004) aboard the Spitzer Space Telescope (| Werner et al. 
2004 ) includes four sensitive, background-limited spec- 
trograph modules (covering 5-38 jum) based on fixed sin- 
gle slits, with no moving parts. The Spitzer telescope 
offers a dedicated spectral mapping mode in which a se- 
lected IRS slit is moved in a raster pattern, stopping to 
obtain one or more exposures at each position, to map 
out extended regions. Since the IRS is a diffraction lim- 
ited spectrograph which delivers a PSF varying in size 
by roughly a factor of 5 over the full wavelength range, 
this mode is very useful for obtaining extended source 
full-coverage spectrophotometry free from strong aper- 
ture bias. Spectral mapping also offers the significant ad- 
vantage of preserving high efficiency for spectroscopic fol- 
lowup of individual point sources, while providing the ad- 
ditional capability of very deep, spatially resolved spec- 
troscopy of large, low surface brightness sources, all in a 
single instrument without moving parts. 

Here we describe an approach to reconstructing full 
spectral cubes from scanned or rastered spectroscopic ob- 
servations using a single slit. In §[2j we outline a custom 
algorithm for cube reconstruction, consider issues of ex- 
tended source flux calibration in §[3l and finally describe 
Cubism, a tool implementing these reconstruction tech- 
niques for IRS spectral maps in §|4j 

2. SPECTRAL CUBE RECONSTRUCTION 

The reconstruction of three-dimensional spectral cubes 
from a suitably obtained set of overlapping input spec- 
tral images (each potentially consisting of multiple spec- 
tral orders) differs fundamentally from image mosaicing. 
The input images carry mixed spatial and spectral in- 
formation, the relative spatial offsets among the input 
spectra cannot be readily computed from the mapped 
data themselves, and ensuring flux conservation requires 
a dedicated approach. Here we outline the general prop- 
erties of an algorithm for spectral cube reconstruction, 
discussing more detail of a specific implementation, along 
with analysis of the constructed cube, in the next section. 

2.1. Resampling Noise 

Modern spectrographs, in particular those which are 
cross-dispersed, rarely form spectral images in which the 
spatial and spectral axes are aligned with detector rows 
and columns. Curved orders which cross many detec- 
tor rows are common. A simple method to deal with 
these curved order data involves first resampling them, 
traditionally with a flux-conserving form of linear inter- 
polation, to create a straightened spectral image. After 
this step, simple summation can be used to extract the 
final spectrum. 

Optimal forms of spectral extraction which employ 
weighted su ms to minimize signal degradation (e.g. 
iHornel IT986D could also be applied to such a straight- 
ened image. In practice, however, an unavoidable diffi- 
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Fig. 1. — An illustration of resampling noise. In the top panel, 
above, an ideal, uniform, curved spectral order, with center indi- 
cated by the gray line, is formed from a flat spectrum source with a 
Gaussian PSF of constant width FWHM=1.5 pixels. It is sampled 
by the pixel elements at middle, and resampled to straighten the 
order using linear interpolation, centered precisely on the known 
order center, below. Extractions of three different pixel widths 
(offset for clarity) are given in the bottom panel, demonstrating 
the purely sample-based noise induced by resampling. 

culty enters with this approach when the spatial profile 
of an unresolved source distri buted a l ong t he slit is not 
well sampled. As discussed by Mukail (|1990l ), resampling 
noise results from the initial straightening step, which 
is retained and amplified by all subsequent processing. 
Figure [T] illustrates this effect using a single, idealized 
curved spectrograph order formed from a flat spectrum 
source with a constant Gaussian point spread function 
(FWHM=1.5 pixels). This continuous profile is then ob- 
served with the detector pixel grid (all assumed to have 
flat and equivalent response), and then straightened by 
resampling using linear interpolation centered precisely 
on the known order center position. In spectra extracted 
from the resampled order data, periodic resampling noise 
of up to 20% is induced in this example, varying in mag- 
nitude with the size of the extraction region. The period 
depends on the rate at which the order center crosses 
row boundaries. Though this noise is mitigated at larger 
extraction sizes, this is no consolation when performing 
further operations on the resampled spectra, such as esti- 
m ating th e row by row slit profile for optimal extraction, 
as iMukal considered, or remapping the spatially resolved 
pixels along the slit to a common sky coordinate system, 
as considered here. 

It should be stressed that resampling noise is fully inde- 
pendent from detector or other instrumental properties, 
or the form of the point spread function, but is rather 
an inherent sampling limitation: information on the dis- 
tribution of flux along the slit has been irrecoverably 
lost. As an example, consider a point spread function 
(PSF) which is much smaller than the size of a pixel. In 
an observed spectral image, you cannot distinguish be- 
tween the spectrum of a point source having been placed 
in the center of a pixel, or near its edge. If interpola- 
tion is performed to shift this spectrum down by, e.g., 
one half pixel, in both cases, two equal interpolated pix- 
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els result. Real spectrographs present varying and im- 
perfectly determined order positions, both of which can 
compound this issue. Better sampling alleviates the ef- 
fect, but it persists at the few percent level even at the 
typical sampling of 2 pixels per FHWM. In diffraction- 
limited spectrographs, the PSF sampling can change ap- 
preciably from one wav elength end to the other. The IRS 
spectrograph modules (|Houck et al.ll2004 ) were designed 
to sample the PSF at the long end of each module's wave- 
length range, and thus are significantly undersampled at 
shorter wavelengths. 

To mitigate the effects of resampling noise, it is vital to 
avoid additional interpolation steps which can introduce 
it. For the construction of spectral cubes, two opportu- 
nities for interpolation can be avoided: initial straight- 
ening of the curved order data, and remapping of the 
rastcrcd set of input spectra onto an output cube grid. 
Although the information lost by an undersampled spec- 
trum cannot readily be recovered without imposing prior 
knowledge on the form of the slit profile, avoiding unnec- 
essary interpolation — especially early in the reduction 
sequence — can minimize the resulting sampling noise. 
In typical spectral reduction algorithms employing multi- 
ple interpolations steps, resampling noise can be reduced 
by ~10 percent using the methods outlined here, strongly 
dependent on the signal to noise of the observations, and 
intrinsic sampling of the spectral instrument. 

2.2. Spatial vs. Spectral Degeneracy 

Creating the three dimensional intensity cube 
S(9, <p, A) from a series of two-dimensional spectra is fun- 
damentally limited by an implicit degeneracy in the in- 
formation present in a spectral image. As illustrated in 
Fig. [21 when spectrally mapping with PSF size compara- 
ble to or smaller than the slit width, an ambiguity exists 
in the interpretation of pixels adjacent in the dispersion 
direction. Given pixel pij in a single two-dimensional 
spectral image, which maps to the spectral cube location 
S(9ki, 4>ki, A m ), the adjacent pixel in the dispersion direc- 
tion Pi+ij could be considered to map either to the same 
wavelength at adjacent spatial positions within the cube, 
S(6k+i,i, <l>k+i,h A m ); or to the same spatial position at 
different wavelengths, S{9ki,4>ku\n+i)- 

In general, both spectral and spatial adjacency are en- 
coded in the slit image. When the PSF is large, very lit- 
tle spatial information remains. When the PSF is small, 
spatial and spectral information are irrecoverably con- 
volved together within the slit. In the latter case, a single 
point source with two adjacent spectral lines just blended 
at the instrument's spectral resolving power will appear 
identical to the spectrum of a source elongated along the 
slit width, emitting a single spectral line. Many image 
slicing integral field spectrograph designs mitigate (or 
hide) this problem by sampling the slit w idth using only 
a single pixel fe.g. iPoglitsch et al.l 12004). though with- 
out methods to scan wavelength at sub-resolution, this 
can adversely impact the reliability of unresolved spectral 
line measurements. Addressing this degeneracy requires 
some assumption, with three approaches possible: 

1. Assume no spatial information is available in the 
slit. Adjacent pixels in the dispersion direction are 
placed into adjacent wavelength planes in the spec- 
tral cube, and all spatial information is provided by 



the mapping observations. 

2. Assume equal spatial and spectral information is 
available in the slit. Adjacent pixels in the slit are 
distributed equally into the spectral cube within 
individual wavelength planes ("alongside"), and at 
adjacent wavelengths ("beneath"). 

3. Assume a varying mix of spatial and spectral in- 
formation, which can be restated as a mixing an- 
gle with respect to planes of constant wavelength, 
varying from 45 degrees (equally spatial + spectral) 
to 90 degrees (fully spectral). 

In this work, we have adopted the first assumption, 
though others may be advantageous in cases of extreme 
undersampling of the slit. 

2.3. Remapping the Spectra 

The reconstruction of a full spectral cube S{0,<f>,\), 
with two spatial (#,</>) and one spectral (A) dimensions 
from a set of two-dimensional spectral images which have 
been obtained by mapping a given region of the sky 
could, in principle, be accomplished in a number of ways. 
Pixels from the spectra could be remapped onto a cube 
grid, and interpolated. Conversely, pixels in an output 
cube could be reverse mapped into the original spectral 
image set, and a weighted average of interpolated val- 
ues drawn. In order to minimize the effects of resam- 
pling noise and to ensure flux conservation, however, a 
method which avoids unnecessary interpolation steps is 
preferable. 

Such a method is available in direct clipped polygon 
resampling, in which the contributing weights of individ- 
ual input pixels are computed by their geometric overlap 
with an input mask and output grid. Polygon clipping 
is a mature technique in the field of computer graphics. 
For the purpos e of spectral resampling, simple clipping- 
algorithms, e.g. ISutherland fc Hodgmanl (|1974h . suffice. 

Since spectral orders are commonly curved, and the 
orientation of the spatial and spectral dimensions can 
vary at each point along the order, a fundamental ele- 
ment in this algorithm is a means of separating input 
pixels by wavelength. For this purpose, an abutting set 
of "pseudo-rectangles" (PRs), each of which outlines the 
region of influence of a single wavelength, can be con- 
structed to cover the curved order. The center of each 
PR is placed at the precise order center, as illustrated 
in Fig. [2) and they may be angled to follow the tilt of 
spectral lines (common in particular in cross-dispersed 
spectrograph designs), which itself may vary along the 
order. In this scheme, an individual input pixel can (and 
often does) contribute to multiple wavelengths. The ex- 
act method of constructing the set of PRs is not critical, 
as long as the central position along the slit is well de- 
termined and matched by the PR centers, there are no 
gaps in coverage, and the pixel sampling of the underly- 
ing detector device is respected (which in general means 
the width of the PR should be one pixel). 

With the set of PRs in hand, the cube reconstruction 
algorithm consists of the following steps: 

1. Clip each pseudo-rectangle against the input pixel 
grid for each of N input spectral images / in a given 
mapping sequence. 
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Step 1 



Step 2 



Step 3 
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Fig. 2. — An illustration the spectral/spatial degeneracy in sampled slits. A monochromatic point source with PSF smaller than the 
slit width of two pixels is mapped in three steps perpendicular to the slit (offset along the slit for clarity). In the resulting spectral 
images, at right, adjacent pixels in the "wavelength" direction can be interpreted as adjacent spatially or spectrally. Single pixel width 
"pseudo-rectangles" , each of which bounds the contributions from a single wavelength, are shown in color. 
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Fig. 3. — An illustration of the two-pass clipping algorithm. From top-left: a single pseudo-rectangle (corresponding to a unique 
wavelength), shown overlain on the input spectral image (red) detector grid, is clipped, rotated and offset into the output cube (blue) pixel 
grid at the appropriate wavelength plane. Multiple overlapping sets of input spectra are accumulated in the same fashion. A contributing 
set of input spectra at a slightly different rotation angle is also shown. 
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2. Rotate and translate the obtained sets of clipped 
partial pixel fragments (grouped by image I and 
wavelength A) into the output cube grid, according 
to position and orientation within the map. 

3. Clip the offset pixel fragments again against the 
output (cube) grid. 

4. Accumulate the overlap area-weighted sum of con- 
tributing pixels in each cube pixel (and, optionally, 
a parallel uncertainty cube). 
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The final area- weighted sum can be expressed: 



(fc) T (k) 



S(&ij , , A m ) 



Sfc=l 12xy ePR„ 

where Ixy 1 is the data value at pixel [x, y] of the kth input 

spectral image, miy is a per-image mask used to disable 
input pixels, and 5 xy ij m is the area of overlap in cube 
pixel [i,j, m] from input pixel [x, j/]. The overlap area 
results from two distinct rounds of polygon clipping: 1) 
clipping the pseudo-rectangle m (corresponding to wave- 
length A TO ) against the input pixel [x,y], and 2) clipping 
the resulting "pixel fragment" (which can be a full pixel 
or a polygonal portion of a pixel), suitably rotated and 
offset, against the output cube pixel A schematic 

of the two-phase clipping approach is shown in Fig. [3] 

The core of the cube reconstruction algorithm is the 
computation of the S pixel overlap weights. Their cal- 
culation depends on accurate positioning of the pseudo- 
rectangles on the input spectrum, including any line tilt 
at that wavelength, and the absolute position of the slit 
center (and thus PR center) as well as the rotation of the 
slit in the output coordinate grid. Errors in positioning 
the PR along the order in the input spectra or spatially 
within the output grid will degrade the achieved spectral 
and/or spatial resolution, as well as increase correlation 
of noise in adjacent pixels. An individual PR may have 
both an effective rotation induced by instrumental line 
tilt in the spectral image, and a true rotation arising from 
differential slit position angle at the time of observation 
with respect to other spectral images in the mapping set. 
Only the latter should be used to rotate the set of pixel 
fragments into the output grid (see Fig.[3|). 

If the spectrograph produces multiple orders from a 
single slit (e.g. using cross-dispersion), the algorithm of 
Eq. [T] can be extended to permit separated groups of 
input pixels in the wavelength overlap region between 
orders to contribute to the same output plane. The 
wavelength sampling of the output cube in the over- 
lap region between orders is chosen to accommodate the 
maximum spectral resolution present there, and the 6 
weighting terms are split linearly between the brack- 
eting output planes, according to their proximity to 
the output wavelength, e.g. 5 — > [fS, (1 — f)S], where 
/ = (A2 — Ao)/(A2 — Ai), for an input wavelength Ao, 
and its bracketing output wavelength pair [Ai, A2]. This 
ensures flux conservation even for wavelengths sampled 
in multiple orders. 

If accurate error estimates are available for individ- 
ual pixels in the input spectral images, they can be pro- 
pogated through Eq. Q] using identical weight factors to 
build a parallel variance cube: 



(i) 



a (8ij,(f>ij, A m ) 



(k) 

where Vxy is the variance per input pixel. This algo- 
rithm does introduce noise correlations between adjacent 
pixels (as would interpolative approaches), but in prac- 
tice these are minimized by employing a high level of 
pixel redundancy in filled spectral maps, with often a 
dozen or more input pixels contributing to a single out- 
put pixel. Such pixel redundancy occurs naturally for 
overlapping spectral orders, and is greatly increased by 
constructing mapping observations in which single posi- 
tions on the sky are sampled by several different parts of 
the slit (along both its width and length). 

The polygon-clipping method discussed here bears 
a superficial sim il arity to the DRIZZLE algorithm of 
iFruchter fc Ilookl (|2002l ), but differs in that it employs 
exact pixel overlap calculations, two levels of pixel clip- 
ping (at the pseudo-rectangle stage, and again into the 
output cube), and does not shrink pixels or pixel fra- 
ments via a "drop factor". Its advantage over related 
methods to remap curved spectral data is primarily in 
avoiding unnecessary interpolation steps, which can in- 
troduce resampling noise early in the reduction process, 
as well as enforcing surface brightness conservation im- 
plicitly. Another advantage is the ability to track indi- 
vidual contributions from spectral pixels in the source 
images to the final cube, which can be used readily to 
reject discrepant pixels. 

2.4. Pixel Rejection and Masking 

The algorithm introduced offers significant statistical 
power for identifying and rejecting transient or persis- 
tent bad pixels. Depending on the sampling redundancy 
at each pixel within the spectral cube, sigma trimming 
methods can be used effectively. Since a given input pixel 
[x, y] will appear many times in a large spectral cube, 
and since each output pixel can have many (dozens or 
more, depending on map layout) contributing input pix- 
els, statistics can be accumulated to flag outlier pixels. 
Various schemes can be employed, but a sigma threshold, 
normalized to an effective measure of the natural spread 
in the data — the median absolute deviation from the 
median — proves effective. To identify persistent deviant 
pixels in all input spectra, it suffices to demand that a 
given input pixel (at a specific [x,y]) is flagged at least 
some fraction of the times it appears in the output (e.g. 
50%). See § 14.51 for a specific implementation of this ap- 
proach. 

3. EXTENDED SOURCE FLUX CALIBRATION 

Accurate flux calibration of extended sources has al- 
ways proved challenging for slit spectrographs, since 
uniform, beam-filling astrophysical reference sources of 
known flux intensity (in MJy/sr, for instance) as a 
function of wavelength are not generally available. In- 
stead, point sources (typically stars) with well-measured 
or modeled flux distributions are used to provide spec- 
trophotometric reference. The flux-calibration of point 
sources is then, ideally at least, trivial — as long as all 
science targets are unresolved, and are placed in exactly 
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the same position within the slit as the reference cal- 
ibration sources, all details of beam acceptance profile 
and diffractive losses, which can vary significantly with 
wavelength and position along the slit, precisely cancel 
out, and the ratio of unbiased spectra between target and 
calibration reference will form an accurate physical flux 
ratio with the model flux. This ideal situation, though 
approachable for point sources, is by definition not possi- 
ble for extended sources, light from which enters the slit 
at all allowable angles, modulated by the source's angu- 
lar distribution within the slit. For this reason, some care 
must be taken to obtain accurate spectrophotometric cal- 
ibration of extended sources, and in particular spectral 
cubes reconstructed from mapping observations of these 
sources. 

There are three main issues with extended source cali- 
bration based on point source spectrophotometric refer- 
ences, which will be discussed in order of their ease of 
treatment. 

3.1. Aperture Correction 

In imaging point source photometry, a crucial issue for 
calibrating extended sources is the aperture over which 
the point source flux is measured, typically stated as a ra- 
dius in pixels. Assuming all unresolved sources are mea- 
sured with the same aperture radius, the derived fluxes 
do not depend on the aperture size. However, since such 
an aperture encircles only a fraction of the total flux of 
a source, an aperture correction must be applied to de- 
rive accurate extended source c alibrations from k nown 
photometric reference stars fe.g. iReach et al.ll2005f ). 

An identical issue occurs in slit spectroscopy, with the 
extraction aperture width serving the role of the pho- 
tometric radius. Small aperture widths are desirable 
for faint sources, since they exclude noisy pixels in the 
source's profile along the slit. Larger aperture widths, 
however, minimize aperture corrections, since they cap- 
ture the bulk of the source flux present in the spectrum. 
This tension between the desire for larger aperture for 
unbiased flux measurement, and smaller apertures for 
high signal-to -noise flux rec overy has lead to PSF fitting 
methods (e.g. lStetsonlll987f ). and their approximate ana- 
log in spectr oscopic applications, optimal extraction (e.g. 
lHorndfl98l . 

Neither PSF fitting nor optimal extraction are appli- 
cable to extended sources, since they rely on the uniform 
distribution of source flux on the detector or in the slit. 
The first step of accurate extended source flux calibration 
is, therefore, forming an estimate of the aperture losses 
as a function of wavelength (the aperture loss correction 
function — ALCF). This function can be applied to the 
normally-extracted spectrum of spectrophotometric ref- 
erence stars to estimate the full, undiminished flux at 
each wavelength in the limit of an infinite aperture. For 
long slits, the form and magnitude of the ALCF can of- 
ten be obtained trivially using bright flux reference stars 
with wide extraction apertures. 

3.2. Slit Loss Correction 

A slit has finite dimensions, and therefore accepts only 
a portion of the flux of a point source imaged onto it. In 
diffraction limited spectrographs, such as the IRS, this 
slit throughput changes appreciably with wavelength, as 
the PSF delivered by the telescope grows in size by 
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Fig. 4. — The slit-loss correction function for the four IRS mod- 
ules. 

diffraction. For calibrating point source spectra, this 
does not present a difficulty; as long as all point sources, 
including the flux references, are well-centered in the slit, 
the slit throughput function divides out. For extended 
source, however, a calibration system tied to reference 
point source observations will necessarily include an im- 
plicit correction for the fraction of the point source flux 
which the slit admits. In the limit of a perfectly uniform, 
slit-filling extended source — which can be conceptual- 
ized as an infinitely dense grid of point sources with van- 
ishing flux — the diffractive losses out of the slit are 
exactly offset by the diffractive gains into the slit from 
sources beyond its geometric boundary. 

In this limit, then, a slit loss correction function 
(SLCF), which is simply the fraction of the flux of a 
well-centered point source admitted by the slit as a func- 
tion of wavelength, must be used to correct extended 
source fluxes. Without this correction, extended source 
fluxes are overestimated. Obviously, the actual magni- 
tude of the true SLCF depends on the observed source 
flux distribution, which can be intermediate between un- 
resolved and fully uniform. This complicates its use, but 
in practice, depending on the size of the slit, this detail 
can often be ignored, and a single correction function 
can be employed for all extended sources. Estimates of 
the SLCF can be obtained by convolving models of the 
PSF (e.g. STinyTim for the Spitzer telescope) with the 
approximate slit dimensions. An example of such a set 
of SLCF's for the four spectrograph modules of Spitzer's 
IRS is shown in Fig. [5] 

3.3. Slit Beam Profile 

Though a slit can be formed with precise geometric 
boundaries, the solid angle it subtends on the sky is 
broadened due to diffraction. Just as for radio dish mea- 
surements, a slit has some beam profile, B(6,<j>), which 
can change with wavelength. Estimat ing this beam 
profile function can be very difficult. Salami ([2000h 
describe a method for de-convolving the spectrograph 
beams aboard the Infrared Space Observatory using a 
series of spectra obtained by rastering a point source 
around the input slit. Such a method obviously relies 
heavily on models of the diffraction of the combined tele- 
scope and spectrograph optical path, but in principle can 
be used to recover B. 

When attempting to test the reliability of a given slit 
loss correction, for instance by comparing synthetic pho- 
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tometry from the spectra to real photometry over a com- 
mon wavelength range (e.g. 24 fixsx Spitzcr photometry) 
there is a strong degeneracy with the effective solid angle 
subtended by the slits, fl c s = J B(6, <fi)d8d(f>, such that 
the SLCF can be scaled up or down, with compensating 
changes to the inferred beam. If the slit is uniform, and 
unless both corrections are separately known, it is useful 
to collapse all information of the beam profile to a single 
number, the effective solid angle per pixel, and derive this 
number by comparing photometry with spectrophotom- 
etry in a sample of extended source measurements. This 
effective solid angle per pixel should be approximately 
equal to the slit width times the pixel scale. 

4. CUBISM 

An implementation of the cube reconstruction algo- 
rithm presented in §[2] for spectral mapping observations 
using Spitzer's IRS spectrograph is available in Cubism, 
the CUbe Builder for IRS Spectral Mapping 7 . Cubism 
is written in the Interactive Data Language (IDL), and 
is designed to combine sets of two-dimensional spectral 
images from IRS mapping observations into single 3D 
spectral cubes, with two spatial and one spectral dimen- 
sion. A variety of cube analysis tools are also available, 
including arbitrary extractions and map creation (e.g. a 
continuum-subtracted line image). 

Cubism, and a manual detailing its use, are available 
from the Spitzer Science Center 8 . 



4.1. Inputs and Calibration 

Cubism uses by default the raw two-dimensional basic 
calibrated data (BCD) spectral images, which have been 
processed by the IRS pipeline (flat-fielded, corrected for 
straylight and certain detector anomalies, etc.). Versions 
of the data without flat-fielding or stray-light correction 
can optionally be used for testing purposes. It applies 
the twin-clipping reprojection algorithm of §[2j using an 
input set of calibration data derived primarily from SSC 
calibration products. A set of pseudo-rectangles for each 
order (see § 12. 3| . which delineate the order position and 
wavelength solution, are constructed directly from the 
order position and line tilt calibrations. The standard 
WAVSAMP calibration files provide a similar set of PRs, 
but do not follow the detector's pixel sampling as re- 
quired to minimize noise correlation, and are therefore 
not employed. Extended source SLCF and ALCF flux 
corrections are applied by default (see §[3]), the latter by 
employing a special form of the standard FL UXCON flux 
calibration input. 

The effective solid angle per pixel (the single free pa- 
rameter in the extended source flux calibration methods 
discussed in §[3j) is calibrated by comparing integrated 
spectro-photometry over ^arcmin 2 areas to corrected 
photometry in MIPS, IRAC, IRS Blue Peak-up, and ISO- 
CAM d ata in bright galaxies from the SINGS Legacy 
survey (jKennicutt et al.ll2003l ). In all cases, these solid 
angle terms agree with the measured slit width within 
10%. 

7 Other tools are available for the pro cessing of single, st aring 
mode IRS observations — e.g. SMART (Higdon et aTJ 20QD and 
SPICE. 

8 http: //ssc . spitzer . caltech. edu/archanaly/contributed/cubism 



4.2. Cube Construction 

Cubism constructs an output cube grid which bounds 
the mapped region, with default pixel size identical to the 
sampling in the relevant spectrograph module (ranging 
from l'/85 to 5'.'08). 

Individual spectra are offset into the output grid ac- 
cording to the requested or reconstructed positions of the 
slit center during the exposure. The latter, based on con- 
temporaneous telescope attitude reconstruction from gy- 
roscopes and star tracker, yield more stable results. The 
limited spatial information available in any individual 
spectral image makes cross-correlation alignment unten- 
able, such that accurate positional accuracy (below the 
level of the pixel sampling) is crucial for this technique. 
Resultant positional errors, compared to fixed astromet- 
ric systems (e.g. 2MASS) are -1". 

Each of the six sub-slits — two each for low-resolution 
modules long-low (LL, 15-37 ^tm) and short-low (SL, 5- 
15 /im), along with two high resolution modules short- 
high (SH, 10-20 ^m) and long-high (LH, 20-38 /mi) - 
are treated independently. A single cube can contain 
data for only a single sub-slit, though for the low- 
resolution modules, it is unimportant how the data were 
targeted (e.g. the second order of SL can be used from 
BCDs in which the first order was targeted, etc.). In the 
latter case, the frame table which defines absolute posi- 
tions of all Spitzer apertures is consulted to calculate the 
position of the offset sub-slit. Consistent with measure- 
ments, spatial distortion along the slits due to projection 
distortion or other optical effects is assumed negligible. 
In each of the high-resolution modules, ten orders are 
formed by cross-dispersion from the single entrance slit. 
The cubes are created by merging the order data in the 
overlap wavelenth region, as described in § 12.31 

The fractional contribution of all input pixels con- 
tributing to a given output cube pixel are computed, and 
are collectively called the cube accounts. These can be 
stored and reread, and permit high level statistical meth- 
ods to be applied to the cube data, for example to exclude 
outliers (see § 14. 5|) . 

For speed, the core clipping functions, which im- 
plement a special-purpose reentrant form of the 
Sutherland-Hodgman polygo n clipping algorithm 
([Sutherland k Hodgmanl 11974ft . is written in C and 
auto-compiled. If auto-compilation fails, an IDL-native 
version will be used as fallback (at a significant speed 
penalty). 

Cubism was extensively tested on simulated spectral 
mapping data to validate the algorithm and ensure flux 
conservation. 

4.3. Uncertainty 

Full error cubes are built alongside the data cubes by 
standard error propagation of Eq. [1] using, for the input 
uncertainty 07 , the BCD-level uncertainty estimates 
produced by the IRS pipeline from deviations of the fit- 
ted ramp slope fits for each pixel. These uncertainties, 
which are statistical ramp uncertainties only, and neglect 
other errors in calibration and later processing steps, are 
used to provide error estimates for extracted spectra and 
constructed line or continuum maps. 

4.4. Background 
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Both to remove the astrophysical foreground and back- 
grounds (primarily zodiacal and cirrus), and to mitigate 
the effects of time varying warm or rogue pixels, it is ad- 
vantageous to subtract near in time background frames 
from each input data record at the 2D level. Suitable 
background frames (which are free from source contami- 
nation in the targeted slit), can be taken from the obser- 
vations directly (e.g. for smaller targets), from dedicated 
offset observations, or, if neither is available, from records 
obtained from the archive at close times (within a few 
days) and ecliptic latitude (\b — bo\ < 10°). If necessary, 
additional correction can be made by subtracting a ID 
spectrum (e.g. extracted from a dark region within the 
cube) — equivalent to removing a constant value from 
each cube plane. 

4.5. Bad Pixels 

The bad pixel masks — m^y of Eq. Q] — are de- 
rived from three sources. Permanently marked non- 
responsive pixels from the PMASK are combined with 
pixels flagged during the exposure for saturation (zero 
or one ramp sample valid), or flat-field difficulties, which 
are reported in the BMASK input frames (automatically 
loaded alongside the BCDs). Additional user bad pix- 
els can be marked either globally (such that pixel [x, y] 
is disabled in all input spectra), or for individual input 
records. Both types of bad pixels can be generated auto- 
matically, by computing outlier statistics in the accounts 
information built up during the cube generation (the full 
record of which input pixels contributed to which output 
pixels, and with what fractional weight). Two param- 
eters, rv, the sigma-trim threshold, and f m im the min- 
imum fraction, control automatic bad pixel detection. 
For each output cube pixel, the contributing pixel resid- 
uals (either with or without subtracted background) are 
tested against an estimate of the contributing pixel de- 
viation a (computed either from a true variance or, for 
small number contributing, the median absolute devia- 
tion from the median value). If a pixel is an outlier, i.e. 
\I xy — median^ygjjm)! > r a o~, it is tentatively marked 
bad. If the same pixel is so marked at least a fraction 
fmin °f the times it appears in teh cube, it is flagged. Sig- 
nificant power at high spatial frequency (point sources, 
edges) results in a large natural range of contributing 
pixel values, such that some care must be taken to ensure 
valid data (e.g. strong spectral lines) are not flagged. 

4.6. Components 

The Cubism interface consists of three main compo- 
nents: the project window, a multi-purpose viewer, and a 
spectrum viewer and map creation interface. Underlying 
the interface components is a scriptable object-oriented 
data container, serialized to disk using IDL's own save 
format (e.g. 'cube.cpj'). All components communicate 
with each other to present a consistent state. For exam- 
ple, when loading a new set of bad pixels in one compo- 
nent, the displayed bad pixels in another are automati- 
cally updated. 

4.6.1. Cubism Project 

All data, calibration inputs, and cube build settings 
are collected and operated on from the Cubism Project 
pane, shown in Fig. [5] From this interface, projects can 
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Fig. 5. — The Cubism Project window, with several disabled and 
selected records. 

be loaded and saved, or recovered from disk. Any num- 
ber of project windows can be opened at a time. Records 
(individual BCDs) can be added and removed, disabled 
from the build, combined into background images in var- 
ious ways, sorted, examined for header information, etc. 
The cube build parameters and cube build itself are con- 
trolled from the project component. The Cubism project 
interface also allows bad pixels and background records 
sets to be loaded, combined, and saved. An individual 
project targets a single slit in a single module. Sepa- 
rate visualization images can be loaded, and feedback is 
provided on the cube build. Various optimizations are 
employed to speed up the cube build process, includ- 
ing keeping track of the specific cube pixels affected by 
changes to bad pixel settings, and re-building only these. 

Using the viewer component, records or stacks of 
records, the assembled cube, or an image loaded for vi- 
sualizing the position of the slit at each step can be dis- 
played. 

4.6.2. CubeView 

Individual records, stacks of records, spectral cubes or 
maps, and visualization images can all be viewed in a sin- 
gle viewer interface: CubeView. CubeView has a num- 
ber of general and special purpose tools that are acti- 
vated as needed. A single project at a time communi- 
cates with one or more instances of the viewer, which 
modifies itself based on the type of data displayed. For 
example, if the cube project sends for view an average 
stack of BCD records, CubeView provides tools for tog- 
gling background-subtraction, editing the slit-end cutoff 
positions of the pseudo-rectangle set, viewing, adding, or 
removing global and record level bad pixels, etc. When 
viewing full spectral cubes, an interface to move through 
the cube planes is provided, and a cube extraction tool is 
enabled, to construct rectangular extractions or matched 
region extractions based on preexisting spectra. When 
viewing cubes, maps, or visualization images, the WCS 
coordinates are reported, whereas with spectral data, 
the wavelength, order, and any pixel flags are shown. 
When viewing visualization images, individual records 
can be selected directly on the overlay (useful for speci- 
fying background records). General use tools, including 
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Fig. 6. — The CubeView viewer, in three different configuration, showing (from left to right), a stack of BCD records with pseudo- 
rectangle set overlayed and bad pixels displayed, a visualization image with the slit positions and selected records shown, and a stack built 
from the cube itself with an extracted region shown. 



histogram scaling, color map modification, line slice plot- 
ting, box statistics, aperture photometry, a pixel table, 
and others, are always available. 

Multiple CubeView instances can be used simultane- 
ously, or a single instance reused to view different types 
of data. Separate projects loaded in the same session 
communicate with their own set of CubeView tools. An 
example of CubeView displaying record spectral data, 
a visualization image, and map created from a spectral 
cube is shown in Fig. [6j 

4.6.3. CubeSpec 

When an extraction is made from a spectral cube being 
viewed, the resulting spectrum is displayed in the Cube- 
Spec tool, shown in Fig. [7] Map images can be assembled 
and directly viewed in the associated CubeView viewer 
by specifying an arbitrary number of foreground and con- 
tinuum wavelength regions over which to average (e.g. a 
line region flanked by two continuum regions), with ei- 
ther wavelength-offset weighted or uniformly weighted 
averages used to set the continuum in each foreground 
plane. Map sets can be saved and restored by name, 
and a custom redshift entered to shift sets (e.g. specified 
at rest frame) into the observed frame of a target. The 
tool also provides simple line fitting capabilities, and the 
ability to create maps using filter curves or other weight- 
ing functions (e.g. a 24 /zm image using the MIPS 24 /im 
filter curve). 

4.7. Outputs 

Cubism outputs data products in standard FITS and 
IPAC Table formats. Although the working storage 
format is IDL SAV file, the full spectral cubes can 



be exported as FITS, with the third spectral dimen- 
sion's wavelengths encoded as a lookup t able, following 
the FI TS spectral coordinate standards of iGreisen et al.l 
(2006). Maps created from the cubes are output as stan- 
dard 2D FITS images, with information encoding the 
foreground and background wavelength range(s) used to 
create them. When uncertainty cubes have been built, 
accompanying uncertainty FITS files are produced along- 
side the FITS data products. 

Spectral extractions are reproduced in the IPAC Ta- 
ble format (an ASCII format with column delineations, 
labels, and units) 9 , and include information describing 
the coordinates of the extraction rectangle. All products 
have associated WCS coordinate systems attached, with 
estimated accuracy of 1". 

5. OPTIMIZED SPECTRAL MAPPING 

As for deep image mosaics, the redundancy of pixel 
sampling is the foremost determinant of the quality of 
a produced spectral cube. Redundancy permits robust 
statistical rejection of bad pixels, minimzes the impact of 
variations in slit throughput or detector response along 
the slit, and reduces noise correlation between neighbor- 
ing output pixels. 

In particular for the IRS, with its time varying rogue 
pixels, sampling a given position on the sky with at least 
4 different detector pixels is crucial for achieving high fi- 
nal cube quality (e.g. using the methods of § 14. 5p . Given 
the spatial and spectral degeneracy (§ 12. 2[) . it is also crit- 
ical to step by one-half the slit width in the dispersion 
direction. Wider step spacing risks gaps or uneven photo- 
metric response in the spectral cube, in particular when 

9 http://irsa.ipac . caltech . edu/applications/DDGEN/Doc/ipac_tbl .html 
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Fig. 7. — The CubeSpec window, displaying an extracted spectrum, and map region comprised of one peak and two continuum regions. 
The line has been fitted. 




Fig. 8. — An example map created from a SL IRS spectral 
map (7-15 £im) of supernova remnant Cas A, in three continuum- 
subtracted lines: [Nell] (red), [Arm] (green), [S iv] (blue). Each 
position in the map has an associated full 5—38 /im spectrum. 

a random and unkown spatial dither pattern is added to 
the requested pointing positions. A detailed set of rec- 
ommendations for planning and executing spectral maps 
is provided at the SSC 10 . 

6. SUMMARY 

In observing regimes with stable astrophysical back- 
ground and image characteristics, integral field spec- 
troscopy of extended sources can be achieved efficiently 
by rastering a single long slit to cover the target area 
with suitable redundancy. This can significantly reduce 

10 http: //ssc . spitzer.caltech.edu/ irs/documents/ specmap_bop 



instrument complexity and size, in particular if a long 
slit mode is already required for efficient followup of iso- 
lated point sources. These savings in instrument com- 
plexity are traded against increased dependence on tele- 
scope pointing accuracy — only telescope attitude sys- 
tems which can track position at angular scales substan- 
tially smaller than the smallest slit width provide effec- 
tive platforms for slit mapping integral field spectroscopy. 

A resampling noise minimizing, two-pass polygon clip- 
ping algorithm for the reconstruction of spectral cubes 
with two spatial and one spectral dimension from spec- 
tral mapping data sets was introduced, and an implemen- 
tation of this algorithm for Spitzer IRS spectral mapping 
observations, Cubism, was described. Cubism has been 
used for individual maps comprised of ^9000 individual 
spectral data frames, for ultra-deep blind spectral map- 
ping surveys, and for large area mosaics covering hun- 
dreds of square arcmin. It makes full use of the unpar- 
alled background-limited performance of the IRS spec- 
trograph to provide sensitive integral field infrared spec- 
trsocopy. An example custom image created from three 
continuum-subtracted lines from a large map of the su- 
pernova remnant Cassiopeia A is shown in Fig. [8j 
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